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We obtain the rigorous WKB expansion to all orders for the radial Kepler problem, using the residue 
calculus in evaluating the WKB quantization condition in terms of a complex contour integral in the 
complexified coordinate plane. The procedure yields the exact energy spectrum of this Schrodinger 
eigenvalue problem and thus resolves the controversies around the so-called " Langer correction" . The 
problem is nontrivial also because there are only a few systems for which all orders of the WKB series 
can be calculated, yielding a convergent series whose sum is equal to the exact result, and thus sheds 
new light to similar and more difficult problems. 

PACS numbers: 03.65.Sq, 03.65. Ge, 03.65.-w 
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In recent years many studies have been devoted 
to systematical investigation of the accuracy of 
semiclassical approximations, which is a very im- 
portant problem, especially in the context of quan- 
tum 

chaos (Casati and Chirikov 1995, Gutzwiller 
1990, Robnik 1998). Among them the WKB ex- 
pansions to all orders for a one dimensional system 
with the potential V{x) = Uq/ cos 2 (ax) was in- 
vestigated by Robnik and Salasnich (1997a) and 
for the angular momentum operator by Robnik 
and Salasnich (1997b) and by Salasnich and Sat- 
tin (1997). An important earlier classic paper is 
by Bender et al (1977). It has been shown for 
the potentials treated in the above papers that 
all WKB orders can be obtained, that the semi- 
classical series for the eigenvalues converges and 
when summed yields the well known exact re- 
sults. In the present paper we apply the same 
WKB method to the radial Kepler problem. We 
solve the problem in a rigorous way, using the 
residue calculus in evaluating the complex contour 
integrals in order to obtain all orders in the H- 
expansion, which confirms the conjecture by Rob- 
nik and Salasnich (1997b) about the form of all 



WKB terms in this problem. When summed, the 
WKB series yields the exact energy spectrum of 
the Kepler problem. 

The problem is very important as it resolves the 
controversies around the so-called "Langer correc- 
tion" (Langer 1937, Gutzwiller 1990), which states 
that in applying the first (leading) order WKB 
method (the torus quantisation of Einstein and 
Maslov) one should replace the exact value of the 
angular momentum L = l(l + l)h by L = (l + ^) 2 h, 
in which case we then get (by brute force) the 
exact result for the Coulomb energy spectrum. 
Of course, as explained by Robnik and Salasnich 
(1997b), this "Langer correction" has no physical 
foundation and indeed is just an ad hoc guess. In- 
deed, they have shown that higher order terms in 
the Kepler (Coulomb) energy eigenvalue problem 
exactly compensate the error done by the "Langer 
correction" , so that after summing up the infinite 
WKB series we get the exact result free of any ad 
hoc assumptions like "Langer correction". Thus, 
the importance of the present paper is to provide a 
rigorous proof of the hypothesis (ansatz) of Robnik 
and Salasnich (1997b) on this problem. We shall 
consider the Schrodinger equation in the complex- 
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ified coordinate space and shall use the standard 
WKB ansatz. Then we solve the quantisation con- 
dition (which is equivalent to the single- valuedness 
condition of the wavefunction in the complexified 
coordinate plane) by using the residue calculus. 
The problem is not trivial, because it is hard to 
guess the correct form of the polynomial terms 
that enter the quantisation condition at each or- 
der, and experience in problems connected to the 
polynomial differential equations (Dolichanin et al 
1998) proved to be crucial. In fact, using this ex- 
perience we went further producing the more gen- 
eral results on WKB method in the Hamiltonian 
systems with one degree of freedom which will be 
presented in a separate paper (Robnik and Ro- 
manovski 1999). 

We consider the Schrodinger equation for the 
radial Kepler problem 



ft 2 d 2 

[- — ^ + V{r)]^{r)=E^{r) 



where 



V(r) 



2r 2 



a 
r 



with L 2 = 1(1 + l)H 2 ,l = 0,1,2,. 
squared angular momentum, and a > 0. 
We can always write the wavefunction as 



(1) 
(2) 

being the 



ip(r) =exp|— cr(r)| 



(3) 



where the phase a(r) is a complex function that 
satisfies the differential equation 

a' 2 {r) + (j\ a"(r) = 2(E - V{r)). (4) 



The WKB expansion for the phase is 

*(»■) = Ely) ^ 



k=0 



Substituting (||) into @ and comparing like pow- 
ers of ft gives the recursion relation 



a' 2 = 2(E-V(r)), 
fe=0 k n—k n— 1 = 0- 



(6) 



The quantization condition is obtained by re- 
quiring the uniqueness of the wavefunction, so for 
the contour integral around a closed contour in 

the complexified coordinate plane r we must 
have 

oo /h ^k 



da = 2_, I — J (f> dak = 2irn r h, (7) 



where n r > 0, an integer number, is the radial 
quantum number and 7 is a complex contour en- 
closing the turning points on the 
real axis. 

The zero-order term is given by 



j dao = 2 J dr^2(E - V{r)) 



V -E 



and the first odd term is 
'fi 



dai = —ttH. 
^From (^) we get 

a' 2 = (L 4 - 6aL 2 r + (3a 2 - 12L 2 £)r 2 + 



8aEr 3 )/(8r(y/2Er 2 + 2ar - L 2 ) 1 



(8) 



(9) 



(10) 



and 



a' 3 = (3r[a 3 r + 2aEr(-5L 2 + AEr 2 )+ 

a 2 (-3L 2 + 4Er 2 ) - 2EL 2 {L 2 + (11) 
8Er 2 )])/(8(V2Er 2 + 2ar - L 2 ) 8 ). 

Computing also a' 4 one finds, correspondingly 
(Robnik and Salasnich 1997b), 



t)/ t - = -^ 



and 



daA = ft, 



64L 3 ' 



(12) 



(13) 



/,From this result the following general formula for 



the ft 2k term 



was conjectured (Robnik and Salasnich 1997b), 
k = 1,2,3,... 

^ ^da 2k = -2vft(^j2- 2k \ l - 2k (14) 
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and so, taking into account that due to the con- 
jecture of Bender at al (1977) 



j> da 2 k+i = 0, 



(15) 



we get 



= = ^r + \) + t(^2-' k X^% (16) 



where A = L/h = \Jl{l + 1), implying the energy 
spectrum for the quantal Kepler problem ([[]), ©, 



E 



-a 



2h\n r + l + l) 2 



(17) 



We will show that formulae ( |14[) and (|15| ) indeed 
hold. To compute the contour integrals (0) 



do-fc 



(18) 



we find the form of the functions a' k . We have from 



the coefficients contains exactly the information 
which we need. It is interesting to mention that 
similar recurrence relations arise in the investiga- 
tion of the famous problem of ordinary differen- 
tial equations, the so-called center-focus problem 
(Dolichanin et al 1998, eq.(15)), and our previous 
experience turned out to be very helpful here. Our 
surmise reads then: 



and 



G 2s 



a 2s+l 



P 



4s- 1 



r(V2Er 2 + 2ar - L 2 ) 6 ^ 1 ' 

rQ^s-i 

(V2£r 2 + 2ar - L 2 f s + 2 ' 



(22) 



(23) 



where s > 0, the degrees of polynomials 
Pis-iiQ&s-i (as functions of r) are not greater 
than 4s — 1, and P± s -l has the form 

Pis-i = ^ (Jj (L is - 6L is ~ 2 asr + ...). (24) 



/ 2Er 2 + 2ar - L 2 

- V 7 2 

V2Er 2 + 2ar - L 2 



(19) 



We see that a' is a two valued function. However 
we can cut the complex plane along the real axis 
between the turning points rx,r% {r\ < r%), which 
are the roots of the equation 



2Er 2 + 2ar - L 2 







(20) 



and in such case we obtain a single valued function 
in the cutting plane. Then from (161) we get 



ar — L z 



2r(y/2Er 2 + 2ar - L 2 ) 2 ' 



(21) 



We now make an "educated guess" concern- 
ing the form of the coefficients a' k . Namely, we 
shall show that the formulae (22)-(24) are crucial 
for our analysis, however it is not obvious from 
the first few expressions for a' k that this shape of 



Here and below we denote by Qk any polynomial 
of degree less than or equal to k. 

by induction 



We prove the formulae (P2[)-(|2l 
on s. According to dl~0|), ( jTT| ) we see that for s = 1 
the statements are true. Let us suppose that they 
are proven for all s < m and consider the case 
s = m. 

To compute a' 2m using the recurrence formula 
@ we note that 



a 2k+l a 2(m-k-l)+l 



r 3 Q<lk-lQ<lrn-4k-5 

(V2£r 2 + 2ar - L 2 ) 6m_1 

_ f 4 Q4m-6 

r(yj2Er 2 + 2ar - L 2 f m ~ l 



(25) 



a 'l°2(TO-l)+l 



0"n 



(26) 



! (ar-L 2 )Q 4 „ 



iV2Er 2 + 2ar - L 2 ) 



6m— 1 
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and 



it 

a 2m-l 



(27) 



rQ 



4m— 5 



\ (V2£r 2 + 2ar - L 2 ) 6m ~ 4 J a' 
(r 2 [(—L 2 + 6ar — 6amr + 
10.Br 2 - 12Emr 2 )Q 4m _ 5 + 
r(-L 2 + 2ar + Er 2 ) x 
QL- 5 ])/(KV2i?r2 + 2ar-L2)^-i). 

Thus we see that the fractions (^)- (P7|) have the 
form (^) and the polynomials in the numerators 
have no influence on the two lowest terms of the 
numerator of cr' 2m . 
Therefore 



°2m 



TO— 1 n-' it' 

y-^ °2fc°2(m-fc) | 



fe=l 



2^o 



(28) 



r 2 Q 



4m— 3 



r(V2Er 2 + 2or - L 2 ) 6m_1 

m— 1 p p I „2 

x - "4fe-l "4(m-fc)-l + r 

~ 2r(V2£r 2 + 2or - L 2 ) 6m_1 



4m— 3 



Note that the degree of the polynomials 
Pik-iPA{m-k)-i is l ess or equal 4m — 2. Hence 
it remains to compute the two lowest terms of the 
polynomial cr' 2m . Observing now that 



Pik-iP^i 



1 



m— k) — 1 



2 2m Ifc H m _ fc 

x(L 4m -6marL 4m ~ 2 + ...) 
we get from (pg|) 



/ 

a 2m 



(29) 



(30) 



L 4r, 



6L 2 mar + HQ4 m _3 
2 2m + 1 r(v / 2£r 2 + 2ar - L 2 ) 6m_1 

m-l /1\ / 1 

xE 

fc=i 

Therefore, taking into account the equality 




(which is just a special case of the Vandermonde's 
convolution, see e.g. (Graham at al 1994)) we con- 



clude that (H) holds. 

It remains to find cr 2m+l . We get 



a \ G 2 m 



(T. 



2 m. 



2a' 



(32) 



((6am + 12Bmr)P 4m _i + 
(L 2 - 2ar - 2Er 2 ) x 



Pi m _ 1 )/(2(V2Er^ + 2ar - L 2 ) 6m+2 ) 
and in the case < k < m we have 



2fc"2(m-fc)+l 



cr, 



(33) 



o 



?"P4fc-lQ4m-4fc-l 

(\/2Er 2 + 2ar ■ 



£2~)6m+2 

Using the induction assumptions (|22| ) and (|2J 
we conclude that ( |32| ) and (|33| ) are of the form 
(23), so each (T 2m+l obeys (|23|). 

We 



Now we can compute the integrals (18). 
are going to use the residue calculus. The inte- 
gration contour in our complex r plane (complex- 
ification of the coordinate space r) is chosen so 
that it encircles the two turning points r\ and r 2 . 
They are connected by a cut so that each a' m is 
single- valued, and has singularities (poles) either 
at infinity or at r = 0. Hence the integration con- 
tour encircles now the poles that lie outside the 
integration contour and, therefore, the calculation 
becomes very easy. Namely, using the residue the- 
orem we get 

2k r 

da 2k = (34) 



2vri 



2k 



- ) (Res a' 2k + Res c 



) fJ 2fcJ 



From ([221), ( g4j ) it is easily seen 



For large r we have according to (| 
vergent series 

i Pdb-1 



(35) 



the con- 
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4fc-l 



a 



(36) 



[l + (6fc-l)( 



Er 2Er 2 



+ ... 



and, therefore, 



Resoo^fc = 0. 
Similarly, we obtain 



j> da 2 k+i = 0, 



(37) 



(38) 



since r = is not a singular point of o-' 2k+1 and 

Resooovjfc+i = 0. (39) 

Hence 



2k 



2vri- 



d<T2fc 



(40) 



Tr 



2k 



-2irh 



1)* ^(_l)fc 2 2fe L 2fc-l ^ 

2 2fc Ui' 



Thus formulae (0), (fBj) and fll6] ) are proven, re- 
sulting 

in the exact energy spectrum (|j~7|). 

In conclusion, we have shown that the residue 
theorem can be successfully applied in comput- 
ing the WKB series of certain potentials in one 
degree of freedom, in particular for the Kepler 
(Coulomb) potential. In this way we have con- 
firmed the conjecture by Robnik and Salasnich 
(1997b) about the form of all the terms in the 
corresponding WKB expansion and thus resolved 
the controversies about the so-called "Langer cor- 
rection" (Langer 1937, Gutzwiller 1990), by ex- 
plaining that by ignoring the "Langer correction" 
and assuming the exact value of the quantal angu- 
lar momentum we indeed get the exact result for 
the energy spectrum after calculating the terms of 
all orders and by summing the WKB series (after 
completing our work the preprint by Hainz and 
Grabert (1999) with similar result (but another 
convergent WKB expansion) was posted). 



It appears interesting to apply this method 
to the investigation of other potentials like, for 
example, the above mentioned potential, V = 
Uq/ cos 2 (a:r). Here the difficulties arise because 
the integrated function has infinitely many singu- 
lar points. 

The most important problem in this context is 
the general formulation of the exact WKB quanti- 
zation condition in terms of contour integrals for 
arbitrary one-dimensional potentials. 
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